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We examine cosmological perturbations in a dynamical theory of inflation in which an Abelian 
gauge field couples directly to the inflaton, breaking conformal invariance. When the coupling 
between the gauge field and the inflaton takes a specific form, inflation becomes anisotropic and 
anisotropy can persist throughout inflation, avoiding Wald's no-hair theorem. After discussing 
scenarios in which anisotropy can persist during inflation, we calculate the dominant effects of a 
small persistent anisotropy on the primordial gravitational wave and curvature perturbation power 
spectra using the "in-in" formalism of perturbation theory. We find that the primordial power 
spectra of cosmological perturbations gain significant direction dependence and that the fractional 
direction dependence of the tensor power spectrum is suppressed in comparison to that of the scalar 
power spectrum. 



I. INTRODUCTION 



Inflation gives a compelling explanation of the flatness, homogeneity, and isotropy of our Universe on large scales. 
It also generically predicts a nearly scale-invariant spectrum of density perturbations, which is consistent with our 
observations of the cosmic microwave background (CMB) and of structure formation. Because of these successes, the 
inflationary paradigm has dominated very early Universe cosmology in recent years. 

In this paper we focus on the prediction of isotropy from inflation. The no-hair theorem of inflation states, roughly 
speaking, that an initially expanding, homogeneous universe with positive cosmological constant, A, and matter 
satisfying the dominant energy condition will become indistinguishable from a universe with de Sitter geometry on a 
time scale of \/3/A pQ. Because of the no-hair theorem, isotropy is generally taken as a prediction of inflation. 

But there could be ways around the no-hair theorem. For example, models with spacelike vector fields that get 
vacuum expectation values can lead to a preferred direction during inflation, evading the no-hair theorem because 
the vector field stress-energy tensor does not satisfy the dominant (or even the weak) energy condition |2j. However, 
such "aether" models have been shown to be unstable [3H5]. 

Recently, another model has been shown to support a persistent anisotropy during inflation [6]. In this model, there 
is a nonminimal coupling between a U(l) gauge field and the inflaton, essentially leading to a time-dependent U(l) 
charge during inflation: 



S = I 0l A Xy 



2 (^)W) " V ^ " ^ V 

Here, the (7(1) field strength, F^, may or may not be the electromagnetic field strength. When the coupling, f(<f>), 
between the inflaton, <fi, and the U(l) field takes a particular form and there exists a nonzero homogeneous U(l) 
seed field, an anisotropy persists throughout inflation even though the space-time is undergoing nearly exponential 
expansion. More specifically, the "electric" field contributes non-negligible extra negative pressure in the direction in 
which it points, which causes space-time to expand more slowly in that direction. 

The model avoids the no-hair theorem by having (1) expansion that is not purely exponential and (2) a coupling 
between the inflaton and other matter. The mechanism for evasion of the no-hair theorem shows up in our results 
in the following ways: (A) all modifications to power spectra associated with the anisotropy go to zero when slow- 
roll parameters vanish and (B) isotropic dynamics is quickly restored if the inflaton-dependent coupling that breaks 
conformal invariance goes to a constant (as is the case at the end of inflation, when the inflaton field relaxes to the 
minimum of its potential) . 

All of the standard energy conditions are satisfied in this model, which means it should not be plagued by stability 
issues as in aether models. The model does, however, suffer from the standard fine-tuning problems of single field 
inflation. Nevertheless, to our knowledge this model could be the first consistent model of inflation that evades the 
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no-hair theorem and includes anisotropy at a significant level. It is therefore interesting to investigate whether the 
model is truly consistent and to investigate its potential astrophysical signatures. 

To that end, in this paper we consider gauge-invariant cosmological perturbations in this anisotropic inflation model. 
We consider and discuss a model generalized from that of [5] , and extend their formula for the relation between the 
anisotropic expansion parameter and the slow-roll parameter to include arbitrary forms of the inflaton potential. We 
also present the dominant effect of the anisotropy on the power spectra of tensor, vector and scalar perturbation 
correlations at the end of inflation. 

Our main conclusions are: 

• The power spectra for gravitational wave and curvature perturbations can develop dramatic direction dependence 
for very small values of the anisotropy parameter 1 if the parameter is nearly constant for a large period of 
inflation. 

• The main cause of direction-dependence of the power spectra is a coupling between the U(l) vector degrees of 
freedom to both tensor and scalar degrees of freedom through the anisotropic background. These interactions 
significantly affect the power spectra of modes after horizon crossing. 

• The ratio of the fractional direction-dependent change in the gravitational wave power spectrum over that of 
the curvature perturbation power spectrum is nearly equal to the tensor-to-scalar ratio. In particular, the 
curvature perturbation power spectrum has much stronger direction dependence than the gravitational wave 
power spectrum. 



For a given scale, the tensor and scalar power in modes with wave vector perpendicular to the preferred direction 
is greater than the power in modes with wave vector parallel to the preferred direction. 2 



• There is no indication that the anisotropic inflation model is unstable, (e.g. There are no ghosts.) This should 
be unsurprising since the stress-energy tensor for matter in the model satisfies the dominant energy condition. 

Many have studied inflationary scenarios with actions similar to ([T]), interpreting F^ v as the standard model 
electromagnetic field strength, in the context of explaining the existence of large-scale magnetic fields in the Universe. 
Initially Parker [7 and then Turner and Widrow [H] showed that magnetic fields produced in an inflationary Universe 
are "uninterestingly small" (i.e., too small to possibly account for the observed large-scale magnetic fields in the 
Universe) unless the conformal invariance of the electromagnetic field is broken. The generation of seed magnetic 
fields starting from the action in ([!]) and a particular f(4>) was considered in [5] and more recently in [10] . Generic 
predictions for magnetic fields in a large class of models, of which the model we consider here is an example, were 
presented by Bamba, et. al. [H] ; the particular realization of the model we consider in this paper is what these authors 
refer to as the "weak coupling case". Magnetogenesis, including the backreaction due to electromagnetic fields, in the 
inflationary scenario we consider here was considered in |12j . For a review of the generation of magnetic fields during 
inflation in a more general context see, for example, [13] . 

More recently, the effect of vector fields during inflation has been studied in the context of their effects on the 
curvature perturbation power spectrum. A "vector curvaton" scenario, in which a vector field with time-varying mass 
and Maxwell-type kinetic coupling term contributes to the curvature power spectrum, was found in [14] to allow 
significant anisotropic contributions to the curvature spectrum and bispectrum if the vector field remains light until 
the end of inflation. A similar massless vector curvaton scenario was considered in |15j and again the possibility of sig- 
nificant anisotropic contributions was found. 3 The anisotropic contribution of vector field perturbations to primordial 
curvature perturbation correlations in various inflationary scenarios was also considered in |16H21j . Perturbations of 
what correspond to our cross polarization gravitational wave degree of freedom were studied in |22j , but in a scenario 
in which a second scalar field, uncoupled to the U(l) field and the scalar field that couples to the U(l) field, causes a 
transition back to isotropic expansion before the end of inflation. 

This paper is organized as follows. In section [nj we introduce the model. In section |ITlj we discuss our philosophy 
and methods for calculating and analyzing primordial perturbation spectra. Finally, in sections IV and [V] we calculate 



the primordial perturbation spectra and briefly discuss stability. We summarize our conclusions in section |VIj 



1 The anisotropy parameter is basically the fractional difference between the rate of expansion in the preferred direction and that of a 
perpendicular direction. 

2 I.e., the parameter g, (see equation l |39[ l), as defined in [2], that characterizes the direction-dependence of the power spectrum due to a 
preferred direction is negative. 

3 Both studies employed the SN formalism in calculating the curvature perturbation power spectra. 
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II. MODEL AND BACKGROUND SOLUTION 

We consider a space-time governed by the following action [6]: 



d 4 x. 



-fJ 



R 

2k* 



(2) 



where g = det(g llv ), R is the Ricci scalar, <f> is the inflaton, and = d^A v — d u A^ is a U(l) gauge field strength. 
For convenience, we'll refer to the C/(l) field as the "electromagnetic" (EM) field, even though it need not be the 
standard model EM field. Here we've defined 



8ttG = l/A/| lanck . 



(3) 



We assume that the background is homogeneous, and that there is a nonzero homogeneous electric field. 4 We 
orient coordinates such that i*y = F m = F vz = and F vx ^q. One could just as easily have chosen to consider a 
homogeneous magnetic field. This choice does not change the form of the background stress tensor, and we expect 
the results of this paper to apply in the magnetic field case as well. However, allowing for both electric and magnetic 
fields of arbitrary relative alignment is beyond the scope of this paper. 

The background space-time is Bianchi I and the metric can be written in the following form by appropriate choice 
of coordinate axes: 5 



ds 2 = a(rj) 2 {—drj 2 + r yij{rf)dx l dx^ , 



where 6 



7x 



- 20M and 7y = for all i + j. 



(4) 



(5) 



Since g is independent of /3, the scale factor, a, completely characterizes the space-time volume. For convenience we 
define a to be the logarithm of the scale factor, so 

a = e a . (6) 

In parametrizing the metric, we've used the conventions of 23J. The solution to the background electromagnetic field 
equation of motion is then [B] , 



PA' 



,-4/80?) 



(7) 



where pa is an integration constant of mass dimension two and a prime indicates a derivative with respect to conformal 
time 77. In these coordinates, Einstein's equations take the form [6] 



o 



12 



P' 2 + 



2 

2 2ai 



p 2 e - 2a -4/3 



2/ 2 (0) 



2a' 2 + n 2 e Za V{(f>) + 



p 2 A K 2 e' 2a -^ 



6/ 2 (0) 



0" 



-2a' p' + 



V 2 A K 2 e- 2a --^ 



3/ 2 (^) 

Given Einstein's equations above, the equation of motion for <j> is redundant. 7 



(8) 
(9) 
(10) 



4 At least we assume that the "electric" field was aligned in our causal patch. We will not consider the effects of regions with differing 
directions of alignment of the electric field. 

5 The form is chosen so that the spatial metric has unit determinant (and therefore scaling or translating 0(rf) does not affect the spatial 
volume element). 

6 An equivalent ansatz would have been: ds 2 = —dt' 2 + an (t) 2 dx 2 + a±(t) 2 (dy 2 + dz 2 ). 

7 Recall that Einstein's equations and the matter field equations are related through the conservation equation, V^T^ 1 = 0, where T\5 is 
the matter stress-energy tensor. 
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It was shown that inflation can occur for suitable initial conditions such that the Universe is initially expanding, 
and that the energy density of the vector field will remain almost constant with respect to the inflaton energy density 
if f{4>) oc e~ 2a [BJ. (Recall that if there's no inflaton-electromagnetic coupling, the ratio of electromagnetic energy 
density to inflaton energy density decays as a -4 .) Let us briefly show how this can occur. 

If expansion is nearly exponential (in cosmic time), then the "slow- roll" parameters, 



dtH 
' H 2 
d 2 H 



v'2 



(11) 

(12) 



are very small compared to one and as usual, H 
expansion is nearly exponential. 

The field equations @, ^ and §L0§, can be cast in the following form 

£' 



2Hd t H 

= £^. 8 Higher derivatives of H must, of course, also be small if 



PA 



2 AI2 



2 2 - 

k p A e 



■4« 



2a 2 f 2 {<j>)a r - 
a 2 K 2 V$) 



a 



where £ 



2e - 6£ 
P'/a'. 



= - 3£ - e£ 



£' 



E^ 
2~d 



:PA 



2e£ - 6E 2 



2e 



-,PA - 6£ 2 



(14) 

(15) 

(16) 
(17) 



The quantities and pa are dimensionless energy densities, normalized by the Hubble scale squared times the Planck 
mass squared. 

In standard single field inflation with an inflaton potential V, for example, one finds from the field equations that 
/- ~ y/f-i so that if expansion is nearly exponential, then the inflaton must be slowly rolling. Taking derivatives of the 

above equations in the isotropic case, one can find expressions for derivatives of V in terms of slow-roll parameters — 

thus yielding requirements of a potential that can give rise to inflation. 
From ( 15 ) and ( 14 ) one finds 



Pa 
Paol' 



djV K<j>' 
kV a' 



2e = 



IPa 
3 a' 



= -4-2 



<W Kfi 
nf a' 



3 - e 
f 2c- 



\PA 



4E = 



3E - eE 



(18) 
(19) 



where 



. £' E" " 



We can glean a fair bit of information from equations (14) - (|T9l without much effort. First, what if expansion 



were purely exponential so that 5 = e = 0? From ( 16 1 we can immediately see that pa and E had better then also be 
zero based simply on the fact that K f 2 , pa, and E 2 are positive. This could be seen as confirmation of the no-hair 
theorem; anisotropy can exist only if expansion is not purely exponential. 9 Similarly, if e is small, then pa and E had 
also better be small. In particular, even in small field models of inflation where typically e <C § <C 1, the anisotropy 
parameters E and pa must be order e or smaller. Second, from ( 18 1 we see that is nearly constant with respect to 



the Hubble parameter if e and E are small. Also from ( 18 ) we see that 



kV a 



-2e + 0(e'/a') + 



(21) 



8 Note that ; 

^=2e(e + S). (13) 
a' 

9 A more direct confirmation of the no-hair theorem comes from supposing <j>' = (and, for simplicity, e << 1) so that V{<f>) functions as 
a cosmological constant. Then from | |16| l and | |15[ 

dlogp A d /V(d>) , . 

at at V 3 



So pa, and thus by \1A\ also e and S, go to zero on the time scale promised by the no-hair theorem. 
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Third, from ( 19 ) , if e and E are small, we see that pa decreases rapidly with respect to the Hubble parameter unless 

r 



< -2 



or equivalently unless 



Now since 



fa' 



djy_ 

kV 



±y/lfay/l -3E/e- 



(22) 



(23) 



(24) 



a ready choice for the coupling function, /, if one wants the energy density of the electromagnetic field (and thus the 
anisotropy) not to decay rapidly with respect to the inflaton energy density is thus 



/(0)=exp|2 CK J 



v 



(25) 



where c is an order one constant. This is the coupling function motivated and examined in [B]. Let's suppose the 
coupling function is of this exact form, so 



Pa 

pACt' 



= -4 - 4c 



= -4 - 4c(2e - 6E 



2e - 4E 



kV a' 
...)(-2e + 0(e'/a') + 



.)- 1 +2e-4E 



(c-l)4-4(3c) - + 



(26) 
(27) 
(28) 



Suppose initially that E -C e. If c < 1 then pa decreases along with E as long as e is small. Anisotropy is wiped out 
(albeit much more slowly than in the case where f(<f>) = 1). If c > 1, then pa initially increases, as does E (see (14)). 



The derivative of the electromagnetic field energy density will thus approach zero, — > 0, and so pa and E will 

(c—l) 

become nearly constant for a time. If E is initially greater than 3e e, then pa and E will initially decrease, (j> will 
climb its potential, and then it will fall back down (slowly) after E has approached a constant [BJ. 

From ( 14 1 one can see that if E is approximately constant then E must be positive. So when the space-time 



undergoes anisotropic expansion in this model (and E is nearly constant) the preferred direction expands more slowly 
than the perpendicular directions. 

When (25 1 holds, we can find an expression for E in terms of the slow-roll parameter during the period in which it 



is nearly constant. Assuming 

0(e) wO(£) c-l>(D(e) 



£<0(e), 



E' 



< 



O(eE), 



/a' < G(e 2 E) 



(29) 



we can set the two different expressions for d^V/V derived from equations ( 18 1 and ( 19 1 equal to each other. Using 
this method we find 



E^ 



1 



1 



Ac 2 



3c 



18c 2 



1 - 2c - 4c 2 
18c 2 



cS- 



assuming c — 1 > 0(e). 



(30) 



The authors of [BJ derived this expression to first order in e for the particular potential V = ^m 2 (j) 2 and argued that 
E generically tracks the slow-roll parameter for general potentials. We find that the expression (30) actually holds 
for any potential V in a slow-roll regime (e, (5< 1). 

As c — > 1, the story is a bit different. For example, if c = 1, looking back to equations (26) - (28) one finds that 
Pa, if it is initially greater in magnitude than C(e 2 ), decreases until it's on the order of e z , and then stays nearly 
constant. From numerical studies it appears that if pa is initially much greater in magnitude than 0(e 2 ), then it will 
rapidly settle to a value much smaller than 0(e 2 ). If the magnitude of pa is initially on the order of e 2 or less, then 
it will stay very nearly constant until the end of inflation. An example with c = 1 is provided in Fig. [l] 
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FIG. 1: Log plot of E and e as a function of e-foldings (Aa = a — an) during inflation. The plot was generated with the 

potential V = \m 2 (j) 2 and coupling function f((f>) = exp - ^ j . The initial conditions were <j>o = 17.5/k, (f>'o = 0, an = —75, 

/?o = and = 0. The constants m and pa were chosen so that initially pa/ p<t, ~ 10~ 6 . Notice that E very quickly settles to 
a value that is somewhat smaller than the square of the slow-roll parameter e. 

The trick of this model is to choose f(4>), given V(4>), such that the electromagnetic field energy density does not 
decay rapidly with respect to the inflaton energy density during inflation. We saw above that a choice guaranteed to 
work is 



f(4>) w exp I 2k 



kV 



d(j> 



(31) 



For example, if V{4>) oc <f> n , then f(cf>) 
we would have 



exp[K 2 2 /n]. What if we were to choose instead, say, /(</>) ~ exp[AK</>]? Then 



Pa 
p A cJ 



= -4-2A— +2e-4E. 



a' 



(32) 



If A is order one, then the anisotropy will rapidly decay. However, if A were large enough in magnitude then the 
anisotropy could persist for a good portion of inflation. 

In our analysis, we will use only the background equations of motion, leaving f(4>) and V{<p) generic. We will then 
be interested in scenarios in which anisotropy can persist over several e-folds — scenarios in which j^r = -2 + 0(e) 
and where pa ~ 9E/2 is approximately constant. We saw that consistency of the background equations and a 
slow-roll scenario dictates that pa must be order e or smaller. We also discussed specific examples of functions, /(</>), 
that can lead to such scenarios (assuming, otherwise, a slow-roll scenario, e, 6 <C 1). In order to calculate primordial 
power spectra, we will use the "in-in" formalism of perturbation theory, assuming 

• e < 1, S < 1 

. p A w 9E/2 < 0(e) 



• PaKPao!) < 0(e). 



III. PERTURBATIONS: SETUP AND STRATEGY 



Our goal is to examine whether the background described in the previous section (slightly generalized from the 
space-time of [6]) is perturbatively stable, and to examine its signature at the level of primordial perturbation spectra. 

We have calculated the quadratic action for dynamical modes in terms of the gauge-invariant variables defined 
in appendix |Xj We calculated the action to quadratic order in perturbations starting with the form of the second 
order Einstein- Hilbert action given in appendix [Bl and a similar expression for the quadratic-order matter action. 
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We worked in Newtonian gauge and used a differential geometry package in Mathematica to massage the quadratic 
action into the (relatively) simple, manifestly gauge-invariant form presented in sections IV and[v| 

Regarding perturbative stability of the background, we find that there are no ghosts (fields with wrong-sign kinetic 
terms), and no other indication of instability at the quadratic level. Here, we take "perturbative stability" to mean 
that dimensionless combinations of fields assumed to be much less than one in the perturbative expansion of the 
action remain small. We find that such small quantities do indeed stay small. 

In the remainder of this section we describe how we set up the calculation and analysis of perturbation spectra; we 
describe the physical scenario, the expression for expectation values in the "in-in" formalism, the definitions for the 



relevant degrees of freedom, and, finally, the current bound on a preferred direction during inflation. In sections IV 
and [V] we calculate power spectra and briefly discuss stability. 



A. Physical scenario 

Perturbations from inflation are usually assumed to be generated in the following way |24) : 

• Quantum mechanical perturbative modes are in their ground state throughout inflation. So the vacuum expec- 
tation value of individual modes is zero, though the variance is generally nonzero. 

• The normalization of the ground states is such that when the modes are well within the horizon, the canoni- 
cally normalized 10 fields, (j), obey a simple harmonic oscillator equation and satisfy the canonical commutation 
relations. 11 

• As modes cross the horizon, their correlations are "frozen in" and translate into classical perturbations that 
lead to, for example, density perturbations that seed the formation of structure in the Universe and lead to 
temperature anisotropies of the cosmic microwave background radiation. 

We shall assume the same, with one complication. We assume the quantity, £ = (3 1 /a! ', which characterizes the 
deviation from isotropy, is nonzero so that expansion of the background space-time is slightly anisotropic, and modes 
that corresponded to scalar, vector, and tensor degrees of freedom in the isotropic background are now coupled. 
(Several scenarios in which this can occur were discussed in section |n ) Because of the coupling of modes, the 
amplitudes of tensor, vector, and scalar perturbations are not separately conserved outside the horizon. As the 
inflaton decays at the end of inflation, the dynamics becomes isotropic again, and tensor, scalar, and vector modes 
decouple. At this point, superhorizon perturbations should be frozen in. We are therefore interested in the correlations 
of perturbations at the end of inflation. Especially if the ?7(1) field in our model were interpreted as the electromagnetic 
field, the details of the reheating process at the end of inflation could also be important in calculating the direction 
dependence of CMB power spectra. In this paper, however, we will only examine the effects of the gauge field on 
curvature and gravitational wave power spectra until just before reheating. 



B. Correlations using "in-in" formalism 



Because in the context of cosmological perturbations as described above we know only the quantum "in" states and 
we're interested in expectation values evaluated at a particular time, we use the "in-in" formalism of perturbation 
theory (see e.g. |25j). We separate our Hamiltonian into a free portion Hq and an interacting portion Hi. The 
interaction-picture (free) fields' evolution is determined by the free Hamiltonian. The expectation value for a general 
operator X at (conformal) time i] can be written as 

(X(r ] )) = (x I (r ] )) + i j v d v ' {[H^x'i^ + iif fd'q' j v d v " ([h^'uh^x 1 ^)}]) + ... (34) 

where the ellipsis denotes terms with more powers of Hj and where X 1 is the interaction-picture operator. 



10 In conformal time, the kinetic term for a canonically normalized field, <j>, in the quadratic action takes the form \<j>' 2 ■ 

11 Specifically, 

[9 V <KV,$),<KVM = -ihS 3 (x-y) (33) 

where r\ is conformal time. 



It should be noted that corrections of quadratic (or higher) order in the interaction Hamiltonian can lead to 
ambiguities when the details of the contour integration are not carefully considered [26] . We will work only to linear 
order in Hj, and therefore we need not worry about such ambiguities. 

C. Decomposition of perturbations 

Since the background space-time is homogeneous, we decompose our perturbations into Fourier modes 

/rl 3 k 
— e^'Sfar,). (35) 

We analyze perturbations about an anisotropic background. Since the background is anisotropic and thus there is no 
50(3) symmetry, perturbations cannot be decomposed into spin-O, spin-1, and spin-2 degrees of freedom and analyzed 
separately. We instead decompose gauge-invariant perturbations according to their transformation properties in the 
isotropic limit. (See appendix [A}) 

There are five dynamical degrees of freedom in our model, corresponding to 

• one scalar degree of freedom, r (spin-0 in isotropic limit), 12 

• two electromagnetic vector degrees of freedom, SA + and SA~ (spin-1 in isotropic limit), 

• and two metric tensor degrees of freedom, E + and E x (spin-2 in isotropic limit). 

In order to analyze the relevant dynamical perturbative degrees of freedom in our scenario, we derived the quadratic 
action in terms of the gauge-invariant variables of appendix [X] Then we eliminated the nondynamical degrees of 
freedom by using constraint equations derived from the action. Finally, we canonically normalized the degrees of 
freedom that correspond to the dynamical "free" fields in the limit as f3' /a' — > 0. Within the "in-in" formalism of 
perturbation theory, we take the interaction-picture fields to be those governed by the dynamics in the /3'/a' = 
limit. 

The quadratic action separates into two uncoupled pieces according to a residual symmetry under parity trans- 
formations. (See appendix |A|) The "odd" sector has two degrees of freedom, E x and SA~ . The "even" sector has 
three degrees of freedom, Er, cL4 + , and r. The fields E + ,E X , and r correspond to fields that are conserved outside 
the horizon during isotropic inflation. Here r is a Mukhanov-Sasaki variable, equal to minus the curvature perturba- 
tion, — C, as defined in, e.g. [37], in a gauge with spatially flat slicing. We will therefore refer to r as the curvature 
perturbation. 



D. Canonically normalized variables 

The canonically normalized fields in each sector ( "even" and "odd" , respectively) are given by 

A+ = f($)SA+ 

A- = f(4>)SA- 

h+ = a(r,)E+/ K and . (36) 
y " ' h x = a(r 1 )E x /n 

f = z(rj) r 

where 

z{rj) = a(?7)~. (37) 
a' 

The fields on the right-hand sides of equations (36 1 are defined in appendix [Aj As mentioned above, in the isotropic 
limit E + , E~ and r are conserved outside the horizon. The other important fact about the fields above is that the 
perturbative expansion of the action is valid when 

E+,E x ,RsA+,RsA-,r«l. (38) 

-Frfx * rjx 



Sec §A22\ in appendix |a| 
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E. Comparison with data 

A formalism for finding signatures of a generic primordial preferred direction in the CMB has been developed [2ll22|. 
In [2] a small direction-dependent contribution to the primordial curvature power spectrum is parametrized by <?* 
where 

P(k)=P (k)(l + g* {n-kf) (39) 

and where h is some preferred direction in the sky. It is postulated that g* will be approximately independent of 
the scale for modes of astrophysical interest and that parity is still conserved. Parity conservation guarantees the 
absence of terms with odd powers of (n • k). Contributions proportional to higher powers of (n ■ k) 2 are assumed to 
be negligible. 

Using this formalism, a nonzero value for was found using 5-year WMAP data at the nine sigma level [35]. The 
central value found for is 0.29 for a preferred direction very close to the ecliptic pole. Since the WMAP scanning 
strategy is tied to the ecliptic plane, this strongly suggests that the nonzero value of g* is due to some systematic 
effect [2H1 HH] • Still, we may reasonably take from the analysis in [35] an upper bound for of 

|S.|<0.3. (40) 

In [30] it is estimated that Planck will be sensitive to values of |<jr*| as small as 0.02. 

Obviously, the gravitational wave power spectrum has not yet been measured, so there is no limit on the analogous 
parameter, <?* gra v, for the gravitational wave power spectrum. 



IV. PERTURBATIONS: ODD SECTOR 



As described in section III C the quadratic action separates into two uncoupled pieces according to a residual 
symmetry under parity transformations. We'll therefore analyze the two "sectors" — which we refer to as "odd" and 
"even" for reasons discussed in appendix [A] — in different sections. We start in this section by analyzing the odd 
sector 13 because it is less complicated than the even sector, having only two coupled degrees of freedom (a tensor and 
a vector degree of freedom) instead of three degrees of freedom as in the even sector. The even sector, which includes 
the curvature perturbation, contains the most interesting physics; analyzing the odd sector is valuable for extracting 
flVgrav and as a warm-up for the analysis of the even sector. 

In this section we present the action for the odd sector to quadratic order in gauge-invariant perturbation variables. 
Then we argue that the form of the action implies that the background is classically stable. Next we diagonalize 
the kinetic term in the action by defining new perturbation variables in terms of which the kinetic term in the 
action is canonically normalized. This diagonalization allows us to identify the fields that should be quantized. The 
Hamiltonian derived from the diagonal form of the action is then separated into a "free" part and an "interacting" 
part, and "in-in" perturbation theory is used to find the autocorrelations (power spectra) and cross correlations of 
the vector and tensor degrees of freedom (see (36)) in terms of the preferred direction and the background quantities 
H and S. The most interesting result in this section is the tensor perturbation power spectrum, given in (80). 

In the odd sector, the action takes the form 



5" 



odd 



"'</ / ^^(-h x *'h x '+-A~*'A-'--h x *h y 
1 (27r) 3 V2 2 2 



2 



4p A a' 2 /3 + \^' 2 (zpa/Z + 6£ 2 - ^Ag£ 2 



; 2 - y + 2Sa 'y + "' 2 (2/5a - 2S + 2A £( S4/3 - S 2 



a'T, + A£c/£ J — iip~h x *A~' + h.c. J J (41) 



where 



'y ^ k^kj — k-^c ^ ~\~ k^c 



2,3 



(42) 



Our odd sector corresponds to the 2d- vector sector analyzed numerically in 1221 . 
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^ (44) 



and /' denotes the derivative of f{4>{r))) with respect to conformal time. Without loss of generality we have set k% = 
and we have taken the preferred direction (the direction along which the background electric field points) to be x . 

By inspection we can see that h x and A" decouple when the wave vector is parallel to the preferred direction (so 
ki =0). This decoupling should be expected due to the enhanced rotational symmetry about the wave vector in this 



A. Preliminary look at stability 

By design, the kinetic terms are canonically normalized. And in the short wavelength limit (k 3> aH), the action 
simplifies to that of two uncoupled harmonic oscillators; there's no indication of instability in the short wavelength 
limit. 

Let's consider the case where k 2 = so the wave vector corresponding to a mode points in the preferred direction. 
In this case, ipL = and Ag = 4. By inspection, one sees that the cross-terms vanish. More explicitly, 



S odd ->fc,-o J dr, J ^(^x^x'+^i-'i-'-^x^x (V 



a 



- \a~*A- [k 2 - y + 2£c/^ + a ,2 {Up A /3 - 2E - £ 2 )) ) . (45) 

When fc 2 - > the action for h x takes the same form as in the isotropic case. Though the effective mass for h x is 
not real for all time (so naively, there's a tachyon), the important point is that h x /a, which we assumed to be much 



less than one in our perturbative expansion of the metric (see (38)), oscillates with decaying amplitude before horizon 
crossing, and then remains constant or decays after horizon crossing. In other words, h x ~ aE x never increases faster 
than a, which is consistent with the perturbative expansion. Similarly, given that2£a'^+a' 2 (14 j S A /3-2I]-I] 2 ) < ty, 

dr) 

p 



the long wavelength solution for A is approximately, A C±f + G%f J j}. Now given that / w a 2 ss H 2 r] 2 



can see that -^-SA ~ (C\ + %a 3 )a 4 (which is decaying) in the long wavelength limit. So clearly the perturbative 

expansion of the action remains valid when fc 2 = 0. 

Now let's consider a wave vector that's antiparallel to the preferred direction, so ki — 0. In this case, ip- = \fpAOt 

and A^ = —2. Then the effective mass squared for h x becomes 

m 2 ff = fc 2 --- a ' 2 (2^ + 9S 2 ). 



Compared to the isotropic case, the effective mass squared for h x receives an additional negative contribution. This 
suggests that h x will grow slightly faster than a outside the horizon. The situation is, of course, complicated by 
the coupling to A~ , but all extra terms in the action when k\ = compared to the terms present when k 2 = are 
small. This suggests that any possible growth of the perturbative fields in this case will be very moderate and does 
not represent an instability. This reasoning will be checked by calculating the power spectra of perturbative fields; 
we can check that the magnitudes of power spectra do not grow rapidly in time. 

The same situation occurs in the even sector; perturbations clearly do not grow when fc 2 = and all extra terms in 
the action when k\ = compared to the terms present when fc 2 = are small. 



B. Diagonalized action 



In general, the canonical quantization of a theory can only proceed once the kinetic interactions have been diag- 
onalized. Usually the diagonalization is accomplished by some constant field redefinition. In our case, we need a 
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time-dependent field redefinition because the "coefficients" in the kinetic portions of the action are not constant. (See 
appendix [Cj) 

The kinetic terms can be diagonalized by performing a time-dependent unitary rotation 



_ f cos^fa) -isin^^(r)f\ fUi 

U 2 



-isint/>g(?7) costp^(r]) 
In terms of the rotated fields, C/j, the odd-sector action takes the form 



_?odd 



drj 



_fk_ 

(27T) 3 



1 (U[V (U[ 



2 \U! 2 



1 [U, 

2 \U 2 



M 



where the Hermitian matrix M is defined 
1 fa" , ,f 



M = r - 



/ 



Ec/' 



f 1 

- 1 E 

fa' 2 



-EAr* — -EAr 2 
2 fc 8 fc 



,^a' 2 (3 + A L 



+ [sin(2-0g)cr 3 - COs(27/)g)cr 2 ] 
+ [cos(2^g)o- 3 + sin^i/^)^] 



JL 

fa' 



2 fc 



2 I / 



a 
a 



— Sa' 



/2 



_l_I S _f E A r 



Ja' ' 2 2 

and where I is the 2x2 identity matrix and we have used the following convention for the Pauli matrices 



o- 2 = 



-i 

1 



and (T3 = 



1 
-1 



(46) 



(47) 



(48) 



(49) 



Physical quantities should not depend on the initial value of ipt- Indeed, we will see that correlations of h x and 
A~ at a time, 77, calculated using the "in-in" formalism of perturbation theory, depend only on the change in after 
horizon crossing. 



C. Correlations using perturbation theory 



In order to calculate correlations, we use the "in-in" formalism of perturbation theory, taking the small parameters 
to be e, S, pA, and E. As discussed at the end of section In] we take 



« 1, 



5 = 



2Hd t H 



« 1, p A « 9E/2 < 0(e), 



p'a 
p A a' 



< 



0(e) 



Given these assumptions and the background field equations (|8j) - ( 10 1 



-2 + 0(e), 
We choose as our free Hamiltonian 



fa' 



fa' 
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2 + 0(e) 



and 



o 



rrodd 

H Q = 



d 3 k 
(2tt) 3 



2 \UL 



U' 



2 \U, 



U 2 



where 



1 V 2 



The interaction-picture fields then obey the following equations: 



drf 



r] z 



U =0. 



(50) 



(51) 



(52) 



(53) 



(54) 
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Each of these fields can be expanded in terms of time-independent creation and annihilation operators as, 



d 3 k 
(2tt) 3 



(27T) 



where the canonically normalized mode functions are 



x (0) (M) = 



^—ikrj 



1 - 



V2fc V~ k v, 

and where the commutation relations of the creation and annihilation operators are 



aU(a^] ={2-K)H i H{k-q) and 



= 0. 



Here, 



(55) 



(56) 



(57) 



(58) 



If we choose /3o = then 7 y (770) = But then if /3 changes during inflation, the coordinates at the end of inflation 
will not be isotropic. On the other hand, if we choose /3q so that /3 = at the end of inflation (when the dynamics 
returns to being isotropic), then the coordinates at the end of inflation will be isotropic. The latter choice is more 
convenient. 

Using the results of the previous section and the form of the matrix M in ( 48 1 , the interaction-picture Hamiltonian 
takes the form 



d 3 k ( 1 (U{^ f 



M^ 1 



ui 



(2tt)3 ^2 \Ui, 
where, 

M« = M - M<°> = A (77, k)l + [sin(2V S )a 3 - cos(2^)a 2 ] f 2 (r), k) + [cos(2^)a 3 + sin(2^)a 2 ] f 3 ( V , k) 
and we have defined 



(59) 



(60) 



h( V ,k) = (7%) -7^%))^ 

/ f 1 



1 (a" f" 



a 



f V 2 



fa' 



1 — -E+ -£A r - -£A, 2 
2 2 k 8 fe 



1 



-p A «' 2 3 + Ah 



= I -4 a' 2 I 1 



fs(ri,k) = 



a' / 

1 (r ^ 

f a 



fa' 



2 fe 



A-i-is- 3 SA r 

fa' 2 2 k 



£A, 2 - -p A a' : 



2 fc 



Our convention for the correlations of the fields will be 



(61) 
(62) 
(63) 

(64) 



where the power spectra are the diagonal entries of the matrix dj. Using (34), the correlations can be written as 



(Ui(p, rj) UM V)) = (P, vWj (q, rj))+i dff{[Hj{rf), U( (p, r,)Uf (a, r?)]) 



More explicitly, the correlations take the form 

Ctj (P, rj) = | X (0) (pr,o , v)\ 2 S tJ + * f drfM$ % r,')I Pno (77', rj) + 

where 

I P (V',V) = f(x ( °W)x (0) *fe?7)) 2 - (x (0) *( P ,r/)x (0) (P^)) 2 



(65) 



(66) 



(67) 
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It is clear from this formula that the zeroth-order power spectra of the fields Ui are isotropic and scale invariant and 
that the cross-correlation vanishes. Here it's convenient to define the function 

i(pv',pn) = ip 2l pW,v) (68) 

where, 

Hx, y) = (—^ - A + — -^2+1) sin ( 2x ~ ^ + (-2 1 2 + 1 - 1 ) cos ( 2x - 2 V)- ( 69 ) 

\2x z y z 2x z xy 2y z 2) \xy xy z x y) 

Solving for the correlations of the variables h x and A~ in terms of the correlations of the rotated variables Ui, we 
find 

P kx (p) = cos 2 ^Cu(fi+sm 2 ^pC 2 2(fi + ^sm(2i>p)(C 12 ($-C2im (70) 

P A _ (p) = sin 2 ^ C n (p)+ cos 2 ^ C 22 ip) - ~ sin(2^ p -) (C 12 (p) ~ C 2 i (p)) (71) 

-CA-K* (P) = C h*A- (P) = c °s 2 ^C 12 (P) + sin 2 ^C 21 (p) + % - sm(2^)(C u (p) - C 22 {p)) (72) 

where we have used the fact that i>_z = ~' l \ } y i - All of the above correlations, and ipp, are functions of time. It is 
understood that these expressions are evaluated at the end of inflation. 
From here on, we will use the short hand notation 

P = Pr,o- ( 73 ) 
Using (66 1 and the expression for M^ x \ the power spectra and correlations are given more explicitly by, 



P%* (P, *?) = \X {0) (P, V)\ 2 +P~ 2 j f fi(v',P) I(PV , ,PV)dv / + (J ' sin (2^rf) - 2^{r,)) f 2 (n> ,p) J(pr/ \pr,)dA 

" cos (2^(r/) - 2^(77)) f 3 (v',p) I(pv',PV)dv'^ j + • • • (74) 
P A - & V) = \X {0) (P, V)\ 2 +P~ 2 { r A(V',P) I{prf,pn)drf - ( I ' sin (2^(t/) - tyM)) f 2 (rf,p) I(pr/,prj)drf 



n 



cos (2^(r/) - 2^(77)) f 3 (r/,p) I(pr,' ,pr,)d V ' j j + . . . (75) 

f, , i/7. ;.,•): ip--\ j cos (2-0p(ry') - 2^(17)) f 2 (v' ,P) Hpv' ,PvW (76) 

sin (2^(r/) - 2^(r?)) f 3 ( V ',p) Iipq' ,pri)dA + ... = -C A _ hx (p, n). (77) 
It's clear from the expression above that the correlations are functions only of the change in the angle tpp. 

D. Discussion 

We are interested primarily in direction- dependent modifications to the power spectra — i.e., modifications of the 
power spectra that depend on the direction of the wave vector, not just its magnitude. Non-direction-dependent effects 
will modify spectral indices, but such effects cannot be disentangled experimentally as due to primordial anisotropy. 
In principle, one could use our method to calculate spectral indices and, for example, relate them to the size of the 
direction-dependent effects. 

The largest direction-dependent contribution comes from the piece involving f 2 . The contribution is given by, 



P~ 2 



(7\m(2^(r/) - WM) MV',P)I(PV',PVW) « -^P- (cos[2^1og(aff/p)] - l) (78) 

assuming is approximately constant throughout inflation, where we've used the fact that ^1 — ~ 3 and the 

relevant integral is calculated in appendix [Dj Modes of astrophysical interest crossed the horizon about sixty e- folds 
before the end of inflation, so for such modes, \og(aH/p) ss 60. 
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When 



f(4>) — exp < 2cn 



kV 



d<t> 



(79) 



for c — 1 C(l) we found that pa ~ ^ e during the anisotropic period of expansion. If the anisotropic period of 
expansion were to last all sixty e-folds before the end of inflation, then we should expect order one direction-dependent 
corrections to the gravitational wave power spectrum for inflationary scenarios in which y/e > ~. Such values of e 
can easily be realized in large-field inflationary models. This analytic result seems to confirm the numerical findings 
in [22]. 

Demanding that the direction-dependent effect on the gravitational wave power spectrum for modes of astrophysical 
interest is less than, say, about 30% would mean that the argument of the cosine function in ( 78 1 is small so that the 



cosine can be expanded in a Taylor series. In this case the power spectrum for h x is approximately, 



(p, v) 



{aHf 
2p 3 



(1 



2^\og(aH/p) 



(aHf 
2p 3 



(1 + 4p A (\og(aH/p)) 2 (l - (h ■ p) 2 )) . 



(80) 



where n is the preferred direction. Thus we may identify 



5* g rav « ~Ap A (log(aH/p)) 2 a -18E(log( aj ff/p)) 2 . 

Note that g* grav is nearly (though not exactly) scale invariant for modes of astrophysical interest. 

Imposing a limit like |<?* gr av| < 0-3 for modes of astrophysical interest corresponds to a limit on pa like 



PA (average after horizon crossing ^5 10 when 



\9* 



< 0.3. 



(81) 



(82) 



V. PERTURBATIONS: EVEN SECTOR 

The even-sector action is much more complicated than that of the odd sector. This sector contains three dynamical 
degrees of freedom that, in the isotropic limit, transform as a scalar, vector and tensor under rotations. This sector 
is further complicated by additional nondynamical scalar variables. 

As in the previous section, we begin in this section by diagonalizing the kinetic part of the quadratic action. This 
process is more complicated for the three dynamical degrees of freedom in this (even) sector than for the two of the 
odd sector, and the smallness of certain background quantities must be exploited; we eventually work in the limit 
^ « e « 1, which is confirmed to be a sensible limit at the end of the calculation. As in the odd-sector calculation, 
we quantize and use "in-in" perturbation theory to calculate power spectra and cross correlations of the scalar, vector, 
and tensor degrees of freedom. The most interesting results in this section are the scalar perturbation power spectrum 
( 111 I and corresponding value for ( 112 1, and also the ratio of the direction-dependent correction to the scalar power 



spectrum over that of the tensor power spectrum (116) 



Instead of presenting the entire quadratic action (as we did in (41 ) for the odd sector), here we present the action 
to lowest order in 6, e, pA, and E. We expand the action assuming that pa, E, and p A /pA are order e or smaller. For 
simplicity, we first present the action to lowest order before elimination of the auxiliary fields $ and VP. (See appendix 
[A] for the definitions of $ and VP.) The action can be written 



^»even 



J d?1 J J^ff [ HtMlH + * f Q H + H t Q t * + $^2*] 



(83) 



where the vectors H and $ are defined by 



H 



A+> 
r' 
h+ 
A+ 
V r / 



(84) 
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and the matrices Mi, M2, and Q are given by 



Mi 



(\ 

















\ 





1 

2 


















1 

2 





















a' 2 - ^ 


2itbLa' 
k 















-2iip'-.a' 






2 '2 
8 K 2^2 PA/ 


V> 









iAV2^' % 





0(e) 



(85) 



M 2 = 



■) - 

oa 



+ (2/3)p A a> 2 ) 



3a z 



9,r 



a 2 k 2 



(2/3)p A a> 2 ) - 
' 2 ) - ^) + ' 




(86) 



Q 



^ 



,0 



1 2 



-3\Z2^'-. 2 + ^S(A r -4) 



V2n Y k 



2 K^Z r fe / 



0( e 3/2) 



(87) 



and V'g is as in (44). Note here the identity 



c/ 2 (A fe --4)/U = -4</f ■ 

Solving the (constraint) equations of motion derived by varying the action with respect to $ and 'J and plugging 
the constraint equations back into the action leads to the action in terms of the three dynamical fields: 



tjreven 



/A?/0[Ht(Mx-Qt^Q)H] 



(89) 



Keep in mind that ip'? is a direction-dependent quantity that varies from zero to plus or minus \fpA, depending 
on the orientation of the wave vector with respect to the preferred direction. The bottom right element of Mi, 
representing (minus) the effective mass for f, is |(~j k 2 ) in the isotropic limit. So if, for example, p A is order 

= 0(e) then we should expect a very dramatic direction-dependent effect on the curvature perturbation power 
spectrum, because the direction-dependent term would be on the same order as the normal, isotropic term (at least in 
the long wavelength limit). In fact, assuming that taking into account the Q^M 2 " 1 Q correction to Mi and properly 
diagonalizing the kinetic term in the action would not weaken the direction-dependent effect on the power spectrum, 
we can get a rough limit on the average value of p A /(K 2 z 2 /a 2 ) during inflation, after horizon crossing. Based on 
the argument of section III E[ we may take a 30% direction-dependent contribution to curvature perturbation power 



spectrum to be an upper limit. Noting that — 



PA 

2e 



2 (2 + 0(e, 5)), the 30% limit translates roughly to 14 
(approximate). 



(90) 



< 10" 2 

average 

Given phenomcnological constraints, it is therefore most interesting to consider scenarios in which pa <C e. Taking 

p A ~ (9/2)E « e, (91) 
by inspection one can see that in the long wavelength limit, 

Q+M^Q = 0(p A /e) (92) 



The first equality can be seen from equations and ||37[l, given that p A must be small compared to K 2 z 2 /a 2 
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and 



Mi = 



(\ 










a' 2 









„/2 fc 2 



\0 *2V2^ »4V2^ 






-i4v^— V'r 

1 iL — fei 

2 ^ 2 



■0(e,p A /e). 



(93) 



We will find, with a careful analysis in the pa *C e limit, that the actual constraint on p^ is much stronger than the 
approximate constraint in ( 90 ) . Thus the pa <C e approximation is valid. 



A. Diagonalizing the action 



Once again, the resulting kinetic terms are not diagonalized and canonical quantization cannot proceed. In the 
Pi < e < 1 limit, the kinetic terms can be diagonalized by performing a time-dependent unitary rotation 



r \ _ ( cos%(?7) -iwa.0Arij\ (U 



A+ 



where 



-ism9 r Xi]) cos6»r(?7) 



I'M = -2V2— & = -2V2- 



Vk 2 



' p A a 



(94) 



(95) 



and where if)'- is the rotation angle in the odd sector, given by (44 1. The rotation of f and A + occurs on a much faster 



timescale than that of h x and A since ip'- = 0{yfpA) and 9L = 0{\J Pa/c)- 
In terms of these rotated fields the even action takes the form 



d>] 



d 3 k 



+ 'u + *' 



-h+'h 



(k 2 - 2a' 2 ) h+fr 



where the Hermitian matrix M is defined 



1 (uiy (u[ 

U 2 



2 \W 



1 (Ui 

2 \U 2 



M 



(96) 



M = (k 2 - 2a' 2 ) I + [sin(2%)CT 3 - cos(26» £ )ct 2 ] I 3-^ j a'' 



(97) 



up to corrections of order e, <5, and p^/e. 15 We've used the same convention for Pauli matrices as in Eq. (49) and, 
again, I is the 2x2 identity matrix. 



B. Correlations using perturbation theory 



The analysis of correlations of dynamical fields in this sector will be very similar to that of the odd sector, up 
to minus signs and replacing tp^ with 9^. It should be noted that the largest direction-dependent corrections to 
correlations in the odd sector are order \/~pa, whereas here we're working to order yj pa/z assuming pa <C e. It 
therefore should be unsurprising that the autocorrelation of the gravitational wave amplitude, h + , has no anisotropic 
contribution at 0(W pa/c)- The same can be said of the cross-correlation between h + and A + . 

Considering now only terms up to order \J pa/c given pa <C e, we choose as our free Hamiltonian, 



rr even 

Hq = 



d 3 k 
(2tt) 3 



-h+'h 



'h+*> 



l(U> 

2\m 



g) + \ (f (vo)k ikj - A 



h+h + * + 



1 fUrV 



2 \U: 



(98) 



15 Recall that, e.g., z" /2z = a' 2 + 0(e, <5) and z'/z = a' + 0(e, 8). 
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where 



The interaction-picture fields then obey the following equations, 



drf 



J V 2 



17/ = 0. 



(99) 



(100) 



As in section |IV| the fields can be expanded into appropriately normalized mode functions and time-independent 
creation and annihilation operators. Dropping terms of order e, pa/c, 8 or higher (including terms with coefficients 
(l^iv) ~ 7 y (Vo))kikj) the interaction-picture Hamiltonian takes the form 



«W&G(3) , " W (3) 



(101) 



where 



M (1) = M — A7 (0) = 3 [sin(26» £ )a 3 - cos(2^)cr 2 ] M c/ : 



(102) 



After computing correlations of the rotated variables using the "in-in" formalism, we can find the correlations of the 
unrotated variables using the equations analogous to equations ( 70 1 - ( 72 ) . 
The correlations are approximately given by 



P?&r,)*\x (0) (p,r,)\ 2 +P- 2 
P A+ &v)*\x (0) (p,r))\ 2 -p- 2 
C fA +(j),v) = ~C A+f (p,r]) « ip- 2 



'sin(2%(r ? - 2%(r?)) S^a'V) I(jn/ )PV )drA (103) 
8in(2%(r ? ') - 2%(t?)) 3^a' 2 (7 7 ') J(p^,p»j)d»/J (104) 



cos(2%(t/) - 2%(77)) 3^ / ja /2 (7 ? ')7(F7 , ,w)d77' 



(105) 



where 7 is defined in ( 69 ) . 

K0 



Assuming and = ^7 are nearly constant during inflation, as in the scenarios we described in section |II| then 

(106) 



a' 

and we may estimate the relevant integral as in appendix ID] Then we see that 



Pf(p,v) 



P A+ (P) V) 



(aHf 
2p 3 

[aHf 
2p 3 



1-2 cos 2-f )log(aH/p)) -1 



1 + 2 cos 2^ log(aH/p) - 1 



C fA+ (p, rj) = -C A+f (p, 77) « sin ( ( 2^ j l.,,(«//..>, 



(107) 
(108) 
(109) 



where -f should be taken as the average value after horizon crossing. 

Now g*, the parameter that characterizes the effect of a preferred direction on the CMB power spectrum, is roughly 
given by 



2 cos 2^ log(a77/p) - 1 



(110) 



The maximal value of for a given wave vector is approximately 2y BA. So even if p^/e is, say, order 10 , 
the argument of the cosine in (110 1 could be significant for modes of astrophysical interest because for such modes 
\og(aH/p) w 60. It's then clear that |<7*| could be order one even for very small values of E and pa- 
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Let's suppose that pA is small enough to satisfy the |<?*| < 0.3 bound of section HIE Then the cosine in (107) can 
be expanded in a Taylor series to give 



Wri)*^ (l-Ui PA (h II ■ i -(h -pi 2 ) 



where n is the preferred direction, and therefore 



.9* 



(in) 



(112) 



Note that g*, is negative, as is g* gr av (see equation (81 1). A negative means that, for a given scale, power is mini- 
mized in the preferred direction. We can understand this general feature in the following way: the pressure contributed 
by the background electric field slows the expansion of the direction along which the electric field points. In other 
words, expansion is slower along the preferred direction. Generically the power in primordial perturbations increases 
in proportion to the Hubble parameter squared; the faster the expansion, the more quickly quantum fluctuations 
are stretched into "classical" perturbations. Since the power of primordial perturbations increases with the Hubble 
parameter, squared, and since in our scenario the space-time is expanding most slowly in the preferred direction, we 
might expect that the power of perturbations with wave vectors parallel to the preferred direction will be smaller than 
the power of perturbations with wave vectors in any other direction. We predict that, generically, models in which a 
preferred direction expands more rapidly /slowly than other directions will lead to positive/negative values of g*. 

The limit \g m \ < 0.3 translates into a limit on the average value of ^ during inflation (after horizon-crossing) for 
modes of astrophysical interest: 



PA 
e 



< 



average after horizon crossing 160(60) 2 ' 

Since pA is assumed to be essentially constant during inflation (as is p$), the limit can be written, 



PA 



average after horizon crossing 



< 10" 



(113) 



(114) 



The measurement of <?* puts a very stringent constraint on the ratio of vector field energy density to the inflaton 
energy density. At the same time, we see that even a very small U(l) gauge field energy density during inflation could 
lead to a significant direction-dependent effect on the curvature perturbation power spectrum. 

Supposing that pa *C e, as we've just seen must be the case in order to comply with observation, the ratio of the 
gravitational wave power spectrum (Pr) to the scalar power spectrum (Ps) is approximately 16,17 



Pt 
Ps 



P 



E+ 



Pp 



Pr 



8^ 



16e 



This fact, in conjuction with (811 and (112), leads to the prediction 



5*grav ^ 1 Pr 
g* ~ 64 /y 



(115) 



(116) 



The direction-dependent effects of a small persistent anisotropy during inflation on the tensor power spectrum are 
suppressed with respect to the direction-dependent effects on the scalar power spectrum by a number of order the 
tensor-to-scalar ratio. This is a consistency condition for the model, given the constraint from observation, pa <C e. 



VI. CONCLUSIONS 



In this paper, we considered gauge-invariant perturbations in a class of models with a persistent background 
anisotropy. After determining the quadratic action in terms of the dynamical fields, we computed the dominant 
direction-dependent effects of the background anisotropy on primordial power spectra. 



In the last equality we used equations and \37\ , given that must be small compared to k?z 2 /a 2 . 

What are identified as tensor perturbations are the amplitudes of the transverse, traceless (TT) part of Sgij/a 2 . We defined &gij.TT / a? = 
2Eij, thus the extra factor of 2 2 . 
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We showed that even a very small persistent anisotropy (with the anisotropy parameter much smaller than the slow- 
roll parameter e) can give rise to a dramatic direction-dependent effect on the primordial power spectra of dynamical 
fields. In an anisotropic background, the coupling between what reduce to the spin-1 and the spin-0 and spin-2 degrees 
of freedom in the isotropic case is extremely important. We showed that such couplings give rise to the dominant 
direction-dependent contributions to the primordial power spectra of tensor and scalar perturbations. 

There has been a fair amount of work on vector fields with time-dependent couplings that are put in by hand, 
assuming exponential expansion. We found that the amount of anisotropy in power spectra are quite sensitive to 
the details of how nonexponential the expansion is, and how long the expansion lasts. Perhaps this sensitivity is 
unsurprising in light of the no-hair theorem. 

We found that for a given scale |fe|, the curvature power, P(k), is minimized when k points along the preferred 
direction. 18 We attribute this feature to the fact that, in the class of models we considered, the preferred direction is 
expanding more slowly than other directions. 

We showed that anisotropic effects are more pronounced in the scalar power spectrum than in the tensor power 
spectra. In fact, we showed that the direction-dependent effects on the tensor power spectrum are suppressed with 
respect to the direction-dependent effects on the scalar power spectrum by a number of order the tensor-to-scalar 
ratio. A priori one might have expected that the tensor power spectra and the scalar power spectrum would develop 
fractional direction dependence of the same magnitude. We find that this is not the case. 

Finally, upon examination of the quadratic action for all dynamical degrees of freedom, we find no indication of 
instabilities in this model. This should not be surprising since the matter stress-energy satisfies the dominant energy 
condition. 

We did not calculate the cross correlation between tensor and scalar perturbations. But one can see from the form 
of the quadratic action 19 that such a nonzero, direction-dependent correlation should exist. The cross-correlation 
effect will be small compared to the direction-dependent effect on the curvature power spectrum, but it could be 
interesting. 
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Appendix A: Parametrization of perturbations 

In the following we use many of the same conventions and notation as in |23j . Since the background space-time is 
homogeneous, we decompose our perturbations into Fourier modes 

S(x\v)= J -^j^Sfaf,). (Al) 

For a given Fourier mode, characterized by the time-independent wave vector fcj, we form an orthonormal basis 
{e 1 , e 2 } for the subspace perpendicular to the wave vector such that 

7 4i e^ = S ab and j^e^kj = 0. (A2) 

Here 7^ is the spatial metric defined in Q. Such an orthonormal basis for the spatial hypersurfaces is uniquely 
defined up to a spatial rotation about the wave vector fc, . To remain properly normalized with the above normalization 
condition, these basis vectors must be time-dependent. 

For definiteness, and without loss of generality, we will take wave vectors to be of the form fcj = (fci,/e2,0). The 
basis vectors can then be written as 



e) = f ) and ej = (0, 0, e?) , (A3) 



- 3 % e^h 



where j^kikj = k 2 



In other words, wc found that g t is negative. 
See equations p5l - p9| ). 
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It turns out that there always exists a choice of basis vectors e\ and e 2 that results in the basis vectors having 
definite sign under what we will call 



k parity: fcj — > —ki 



(A4) 



Our basis (A3 1 is such that under k parity, ef —> (— I) a e". Such a choice of basis is now unique up to discrete spatial 
rotations around the fc; axis by multiples of n/2. 

We parametrize the most general perturbations about the background Bianchi I metric Q in the standard way, 



Following 



ds 2 = ~a{ri) 2 [(1 + 2A)drf + 2B { dx l dt + (7^(77) + /i^efeW] 
Bi = diB + Bi 

hij = 2C ( 7ij + ^) + ^djE + 2%^ + 2E tj 



where cry = and H = a' and also, 



-■ ■'),!{, 0. -■■<),!■:, 0. -■•<), I:,, I) and ^ E i: 



0. 



(A5) 

(A6) 
(A7) 

(A8) 



We parametrize perturbations of the inflaton field and the electromagnetic field by 8<f) and SF^ Ul respectively. 
One can show that the following are U(l) gauge and diffeomorphism invariant variables, 



$0) = A 



1 



#(fc) 



-C 



0(77) 
a'(ri) 



a(ri) 



a I^B - 
B 



(k 2 E)' 
k 2 
(k 2 E)^ 
k 2 



E, 



v > 



X{k) = 5<f> + <f>'(r)) 



B 



(k 2 E)' 



$( j (k) = 5F ij +2F v[ iik j 



B 



(k 2 E)' 



*f (fc) = 5F vi - ^ k F ri3 ih(ik k E + E, 



J- n 



B 



(k 2 E)' 



(A9) 

(A10) 

(All) 
(A12) 

(A13) 
(A14) 
(A15) 



The perturbation in the gauge field can be decomposed along directions transverse and parallel to the spatial wave 
vector: 



SAi = (iSA^'+\k,v))el + (5A^-\k, v ))e 2 + (iSA^tj)))^ 



(A16) 



where the amplitudes SA^'^'^k, 77) are U(l) gauge invariant. 20 In Aq = E = B = Bi — gauge the electromagnetic 
gauge fields 8A^- L ' ± \k,ri) are simply related to the gauge-invariant magnetic and electric field perturbations. In 
particular we may define 



M + (M) = ^ ^ ' J 



and 



5A-(k,r,) = - 



k 2 



(A17) 



where j^kikj = k 2 and where spatial indices are understood to be raised and lowered with the spatial metric, 7^. 
The dynamical, gauge-invariant dynamical electromagnetic variables are 5A^(k, 77) as defined above and are equal to 
<L4.(- L ' ± ) (A;, 77) as defined in (A16l in A = E = B = Bi — gauge (a modified Newtonian gauge). 



The factors of i accompanying some perturbations is to ensure that the relation 8*(k,rj) = <5(— k,r]) holds for all Fourier amplitudes. 
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The tensor perturbations, Eij are gauge-invariant by construction. We will further decompose the tensor pertur- 
bations by constructing the two independent symmetric traceless tensors that are transverse to the wave vector 
We again follow and define these tensors as 



E lJ =E+e++ i E x ef p (A18) 

(A19) 



11 9 9 

p.p. — p.p. 



V2 



(A20) 



We have chosen this normalization since, 



7 < V'4^=* AA '- ( A21 ) 

Because we have chosen a basis with the property that, under k-parity, e° — > (— l) a ef these tensors have k-parity 

transformations ef, — > +e% and ef„ — > — ef.-. 

v *j *j £ j 

We will take the Mukhanov-Sasaki scalar variable (which is conserved outside the horizon in the isotropic limit) to 



be 

a' 

r= 17 x + *- (A22) 
9 

In a gauge with spatially flat slicing, this variable corresponds to minus the curvature perturbation, — f, as defined, 
<■!!■■ in EZ). 

Some of the variables listed are not dynamical and must be removed from the action using constraint equations. 
There are a total of five dynamical variables in the theory. In the isotropic limit, these variables correspond to 
two electromagnetic perturbations, two tensor perturbations and one scalar perturbation. Furthermore, the action 
separates into uncoupled parts according to the transformation of fields under k parity: a piece including E + , 8A + 
and r and one including E x and 6 A~ . 

Appendix B: Quadratic action and Einstein's equations 

Given a metric g p „ = <? M „ -I- Sg^, the Einstcin-Hilbert action to quadratic order in Sg^ can be written as 
S (2) Seh = J ' d 4 x^{^9^(y a Sg^)^Sg av )- _L^(V«^)(V^^) 

+ ^5^^ CT (V Q %,)(V Q <5 ffpCT )-^5^^(V«^ p )(V Q ^ CT ) (Bl) 
+7^R^9 pa (Sg w )(6g^) - ^R» v g»°(8g^)(5g p(T ) + ^R{g^5g^) 2 - ^Rg" u g p °{8g MP 6g vl ,)} 

after dropping boundary terms. In the above equation, the covariant derivatives (V) are compatible with the back- 
ground metric 

Vaff^ = 0. (B2) 

We used this form of the action and our parameterization to compute Einstein's equations. In particular, the 
first-order change in the components Einstein tensor can be written in the following way (in Newtonian gauge, where 
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E = B = Bi = 0), 



a 2 <5G^ 



a 2 <5G'' = 



-2A# + 67^' 



ft 



■$0,* - + m'aj + (6ft 2 - a 2 )$ - -(a 2 



ft 



,2 w i - 



2^^ + 0*^ + 3^^* 



ft J 



a 2 5G* = <5] [2*" + (2ft 2 + 4ft') $ + A($ - *) + 2ft<I>' + 4ft*'] - d l <9 3 ($ - $) - 2°^-d j) d k ty 



-ft 



ft 2 



ft'\ T A* . 

— - \I> H $' 

ft 2 / ft 



CT 2 ($ + (*/W)') 



ft 



(£})" - A£* + 2ft(£j)' - a l k {Ef )'5) + 8){o$ d k &) - 2U 1 lk d {k <^ 



i) 



ft' *' 

2^-2--2(4> + *) 



ay/ 

2 — * -4ft$ 
H 



+ 2ft [aidjE k - o k d k E l ] + [((T i k yd j E k ~ (a k )'d k E l ] , 
where ' denotes derivatives with respect to conformal time and 

a' 1 , 



ft 



(B3) 
(B4) 

(B5) 



(B6) 



In these equations, spatial indices are raised and lowered with 7^ . 

Our expressions (B3) - (B5) match those of [23] up to factors of the anisotropic stress, <rj + 2ft<rj, which in [23] 
was set to zero. Note that the Einstein tensor is gauge covariant rather than gauge invariant. 



Appendix C: Diagonalizing a kinetic term 

Suppose a kinetic term takes the form 

K = \x v X' + X'U.X + .XM'.X' 



(CI) 



where X is a vector of fields and M is a time- dependent matrix. Diagonalizing the kinetic term requires a change of 
variables 



X — > VY, 

where V is a time-dependent unitary matrix, such that 

K — ► -Y^'Y' + total derivative + Y^QY 

where Q is some Hermitian matrix. We can calculate directly that 

K = \y V Y' + (v v V + 7+ (Aft - M)F)Y' + total derivative + Y^QY. 

The kinetic term is diagonalized by a unitary matrix V that satisfies 

y+'V = -\rf(Mi -M)V or equivalents VV V = M - M+. 



(C2) 
(C3) 

(C4) 
(C5) 



If M were a time-independent matrix, then the kinetic term would be diagonalized by a constant unitary matrix 
V such that 



0(M- Aft)V = D 



(C6) 



where D is a constant diagonal matrix. 
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Appendix D: Estimates of integrals 



In order to get a quantitative estimate of the effect of the anisotropic background on power spectra, we must 
estimate the integrals in (74) - (77). We may take pa, S, and the slow-roll parameters to be nearly constant. Then 
the relevant integrals are 



p- 2 [ sin (2V>j?(?/) - 2^%(ry)) a'(r)') 2 I( PV ' ,pr,)d V ' , p~ 2 [ cos (2^(V) - 2^)) a'(^') 2 Hpv' ,PVW (Dl) 



"(^(V) - e 2n0 ^)I(pr]',Pv)dv' and [ a'(r)'fl{pn',pr))dr)' 



(D2) 



where I(x,y) was defined in (69) as 



s / i i 2 i i\ , / i i i r, 

HP, V) = ~^2 + ^ + O Sin 2X - 2 V)+[-2 2+ C0S ( 2X - ^ ■ 

\2x z y z 2x z xy 2y z 2 / \x z y xy z x y ' 

During slow-roll inflation, 



a'{rj) = e a{v) H( v ) « -i 



(D3) 



(D4) 



k 2 e~P° — 

M 7 ?') ~ ^) ~ ("W) - at(r))) — V/3a 



(D5) 



( e 2«/9(V) _ e 2«/3(„ )) „ 2nS ( a (,/) _ q,^)) . 



Let us define a new variable z by 2 



prj = e 



(D6) 
(D7) 



From (D4) it's clear that 



clH 

e 2 w and so z w log(H/p) + a. 



(D8) 



We may thus rewrite the integrals (Dl ) and (D2) in terms of the variable z: 



Is^p- 1 I Bin (2^(z I(-e- 2 ,-e- 2 *)e 2 dz, J c = p" 1 j * cos (^(z - z») ) /(-e^ 2 , -e~ 2 * )e 2 dz 



A^p- 1 / *(z-z )i : (-e- 2 ,-e- 2 *)e- 2 dz, ^P -1 / * /(-e^ 2 , -e~ 2 * )e 2 dz 



(D9) 
(D10) 



where z* is the value of z at the end of inflation and 



V^r = p 2 e-P° 
a' p 



Pa- 



(Dll) 



The function 



7(-e- 2 ,-e- 2 -)e 2 



(D12) 



This is just a convenient dimensionless variable and is not equal to a<j>' /a' as in (|37| 
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FIG. 2: The function e z I(— e z , — e~ z *) on a linear scale. The axes cross at the point {0,0}. For < z < z, the function is 
well approximated by — |e 2z *. The frequency of oscillation for z < does not vary much as z„ increases — only the amplitude 
changes. The plot above was generated using z* — 15. 

oscillates rapidly with growing amplitude for z < 0. See Fig. [2j For z > and values of z* on the order of tens, the 
function is well approximated by a constant 



< Z < Z*. 



(D13) 



The constant can be found by expanding the function about z* = oo and then about z = oo. 

The contribution of terms that go like I\ will be subdominant compared to contributions from terms proportional 
to the other integrals 22 , so we will not bother to calculate I\. Since the dominant contribution to the other integrals 
will occur when z > (which corresponds to after horizon crossing) we may approximate the integrals by 



„2z, 



p-^* sin ( 2 f (*-«.))* 



-e "p 



\ a' 



2^ 



(D14) 



./o 



^2^(z-z,)^)dz 



sin — - z a 



(D15) 



/ 2 =^- 1 / 7(-e 



)e 2 <iz sa --e 2z *p 
o 



(D16) 



Modes of astrophysical interest crossed the horizon about 60 e-folds — plus or minus a few — before the end of inflation. 
Such modes of astrophysical interest therefore correspond to z* 60. 



The contribution from 1\ can be important if inflation lasts a very long time — on the order of 10 e-folds. 
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